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Salt-induced aggregation of 20 nm colloidal silica is followed by light transmission, which shows 
an a kinetic form exp[— (t/to) a ], where a = 2.6 and to is an empirical time constant which reflects 
the colloidal stability. We found a power law dependence of to on ionic strength, which can be ex- 
plained by the classical DLVO theory. The neutral polymers polyethylene glycol (PEG) accelerate 
the aggregation rate, and those with higher molecular weight are more effective in inducing the ag- 
gregation with similar stretched exponential form of kinetics. Current theories of polymer-mediated 
interactions provide a reasonable interpretation of the effect of PEG. The stretched exponential 
kinetics of the light transmission is found to be consistent with a cluster-size dynamic scaling model 
of aggregation. 



PACS numbers: 64.70.Nd, 64.75.- 



I. INTRODUCTION 
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If a particle in solution is small enough, the gravita- 
tional potential energy over a macroscopic length is less 
than the thermal energy. Such colloidal particles dis- 
tribute themselves evenly in a vial by thermal diffusion. 
If the colloid were to have no surface charge, the parti- 
cles would collide, stick together, and eventually deposit 
at the bottom of the vial. Adequately charged colloidal 
particles repel each other and stay in suspension essen- 
tially forever. This property of colloidal stability is of 
immense practical relevance to real world materials such 
as food, adhesives, cosmetics, inks, and paints. The silica 
nanoparticle is of longstanding technological significance 
[l| and in a sense is the ancestor of the modern nanoparti- 
cle. Today it has uses ranging from chemical- mechanical 
polishing of silicon wafers to serving as a DNA carrier in 
non- viral transfection Q • 

Colloidal silica aggregates with the addition of salt, 
which lowers the particle surface potential and Debye 
screening length. The presence of a non-adsorbing poly- 
mer can also induce aggregation. This is understood as a 
"depletion force", "entropic force", or "molecular crowd- 
ing" generated by a contest for space between the colloid 
and the free polymer. Asakura and Oosawa (AO) [!, 0] 
were the first to interpret this simple physical concept. 
Further theoretical progress has been made over the past 
decades, including those using the methods of integral 
equation theory [5| and the scaling/renormalization ap- 
proach [1, 0, 1, 1 . 

Here we study the destabilization of colloidal silica by 
both monovalent salt and the polymer PEG. We should 



mention here that PEG is also of ascending importance 
in the science of biomaterials. In enzymology, it has 
been used to study hydration effects, steric hindrance, 
and molecular crowding. When covalently bound to a 
surface it serves as a biocompatible passivation layer. 

Charged colloidal particles in saline solution are typ- 
ically modeled as bodies with repulsive (stabilizing) 
Coulomb interactions and attractive short-range disper- 
sion forces, an approximation known as the classical Der- 
iaguin, Landau, Verwey and Overbeek (DLVO) theory 
flfjl . [Fy]. The electrostatic potential is modeled by the 
Poisson-Boltzmann (PB) equation, which describes the 
relation between charge density and electrical potential, 
under the assumption that the saline ionic charges de- 
viate from their bulk concentration according to Boltz- 
mann's law. For moderately charged spherical parti- 
cles, the external electrostatic potential is of the form 
4>o exp(— rn)/r, where 1/k is the Debye screening length. 
In a monovalent salt solution, k 2 — e 2 n/ ee^kBT where 
eeo is the solution dielectric constant and n the num- 
ber concentration of salt. In aqueous 1 M monovalent 
salt solution, for example, 1/k = 0.3 nm. Increasing 
ionic strength also lowers the particle surface potential 
<f>o- These two factors combine to decrease the electro- 
static interaction potential U e i e and hence lower the sta- 
bility of a colloidal suspension. 

Attractive dispersion forces are the result of induced 
dipole-dipole interactions; these are responsible for the 
irreversible binding in aggregation. For two volume ele- 
ments of material the interaction force is proportional to 
A/r 6 where A is the Hammaker constant. The resulting 
potential is also called the van der Waals potential U v dw- 
The total interaction is then U e i e + U vc iw- 



2 



With increased ionic strength the colloidal particles 
undergo irreversible formation of larger and larger clus- 
ters, aprocess known as aggregation. It has been shown 
[l2l [l3f that the cluster kinetics, size distribution as 
a function of time, and cluster fractal dimension fall 
into two universal classes: diffusion limited aggregation 
(DLA) and reaction limited aggregation (RLA). In RLA 
single particles and clusters collide many times before 
joining together irreversibly. In the DLA limit, parti- 
cles and clusters aggregate the first time they collide, so 
particle diffusion limits the rate of aggregation. The clus- 
ters generated in these two regimes have different qualita- 
tive appearance; the fractal dimension of DLA particles 
is d f ~ 1.7 and for RLA d f ~ 2.0 0, GS El- These 
two classes also exhibit different time evolution scaling 
laws in cluster-size distribution [12, [l3| . Cluster growth 
had been experimentally observed for quite some time 
[l4| , but did not have the benefit of a modern conceptual 
framework. 

As the clusters have fractal dimension less than 3, it 
is possible for a small concentration of particles to span 
the sample volume. The cluster radius of gyration R g = 
bn 1 ^ d f where b is the monomer radius and n is the number 
of particles in the cluster. As the bulk monomer density 
approaches the cluster particle density n/B? g , the clusters 
in effect span the sample volume. 

There are similarities and differences between aggrega- 
tion and a second-order phase transition. As the process 
of aggregation proceeds, the particle size becomes larger 
and larger. In a second-order phase transition, the cor- 
relation length (analogous to cluster size) diverges, but 
the correlations are dynamic, in that particles are free to 
leave and rejoin large clusters, whereas in aggregation the 
particles irreversibly join a cluster. A salient feature of 
increasing cluster size is increasing light scattering. The 
scattered light I{k) ~ n 2 k~ d f for k > R^ 1 [J] where k is 
the scattering wave vector related to the scattering an- 
gle. The total scattered light (I(k) integrated over k > 0) 
determines the turbidity of the sample, and it has been 
used to characterize the correlation length at the critical 
region in colloid-polymer phase separation [l5j . 

The surface of silica in basic aqueous solution is nega- 
tively charged. Silica (SiO^) reacts with water to create 
silanol (SiOH) surface groups [lq . At a high pH, protons 
are pulled from the surface, leaving negatively charged 
silane groups (SiO~), thereby lending colloidal stability. 

A table-top example of the above concepts is provided 
by a 4 % w/w solution of 20 nm colloidal silica at pH 10. 
As the monovalent salt KC1 is increased past 200 mM, 
the colloidal particle surface potential and electrostatic 
screening length decrease, and the silica particles begin to 
combine into larger and larger clusters. As a result, the 
turbidity (light scattering) progressively increases. In the 
presence of PEG, the turbidity occurs more readily. We 
elaborate on this simple demonstration to further our un- 
derstanding of colloidal stability and polymer-mediated 
interactions. 

The rest of this paper is as follows. We first continue 



with a more detailed account of the necessary theoreti- 
cal background. After a description of the experimental 
materials and methods, the data is presented and inter- 
preted in terms of DLVO theory and the two-particle 
Smoluchowski equation, the depletion force, and the dy- 
namics of cluster-size distribution in the RLA limit. 



II. THEORETICAL BACKGROUND 

A. The surface charge of colloidal silica 

Acid titration data of colloidal silica provides some sur- 
prising details about the nature of its surface [13, ITsj l- 
The results can be qualitatively explained as follows: At 
increased pH, the chemical potential of the solution pro- 
ton gas is lower, and the silica surface protons are able 
to enter the solution, leaving behind a more negatively 
charged surface. With increasing ionic strength, the elec- 
trostatic potential for a given surface charge is less, so the 
particle can assume more charge for a given pH. The re- 
lationship between surface charge, surface potential, and 
ionic strength is explained with the PB equation. How- 
ever, it does not explain the titration data. By proposing 
a surface capacitance, the Stern model is successful in 
modeling the titration data. The details of this calcula- 
tion are in Appendix [A"l In Fig. [1] we present the results 
of this calculation for a solution of 4 % w/w 20 nm silica 
with 140 mM ammonia buffer. It should be emphasized 
that this family of curves accurately models the exper- 
imental titration data [l?], EH for colloidal silica. The 
success of the Stern model in explaining this titration 
data lends confidence that we know reasonable values for 
the surface charge density. 



B. Two- particle interactions 

The two-particle interaction is modeled with a electro- 
static repulsion and a short-range dispersion attraction, 
and the polymer-mediated interaction is treated as an 
AO depletion attraction. 

We use estimates of the two-particle electrostatic inter- 
action U e i e based on linear superposition approximations 
[Sll: 

U e i e (h) = 327ree (fcT/e) 2 a[tanh(e^ rf )] 2 exp(-hn) (1) 

where a is the particle radius, r is the particle center- 
to-center separation and h is the surface-to-surface sep- 
aration r — 2a. cf>d is the electrostatic potential at the 
particle-solution interface. An essentially identical in- 
teraction is obtained following the work of Behrens and 
Grier [21]. In the Appendix it is shown how to obtain 
the surface charge density a and hence (fid from given pH 
and ionic strength. 
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For two spherical particles, the van der Waals attrac- 



tion potential is U v 



dW 
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where A is the Hammaker constant [10]. This formula 
results from the integration of the A/r e force described 
in the Introduction. 

For 17jj we use the results of the PRISM theory ac- 
cording to Fuchs and Schweitzer [f| 



U D = k B T 



27 c a 



[1 + (5/9)Jf + (X/3) 



(2) 



where X = h/£o. £o is the polymer mesh length or equiv- 
alently the polymer density-density correlation Q . In the 
dilute limit it is the Gaussian radius £o = Rg/V%- We use 
radius of gyration values R g = 11 A nm for 35 kD PEG 



and R a 



3.1 nm for 4 kD PEG. c* marks the semidilute 



polymer concentration where the molecules begin to over 
lap, i.e. c> 1/(4ti\R3/3). We use c* values of 0.92 % w/w 
and 5.25 % w/w respectively. 



C. Aggregation and the Smoluchowski equation 

The aggregation of colloidal particles is a stochastic 
process. For a spherically symmetric configuration the 
Smoluchowski equation [lOj, |20j provides a continuum de- 
scription of net flux of two particles: 



1 Or 2D Or 



(3) 



where r is the center-to-center particle separation, n(r) 
is the particle density, D is the diffusion constant, and 



potential. This equation is of the form y' + by = f(x) 
and can be solved by using the integrating factor e bx . The 
boundary conditions are n(2a) — and n(oo) = no- That 
is, the particles fall into a sink when they touch, and at 
a large separation the particle density is the bulk value. 
Assuming the quasi-equilibrium condition of constant </, 
the specific solution is (loj 



J 



AttDhq 



JZ + s^[U(r)/k B T]dr/r2 



(4) 



Here 5 is an arbitrary cutoff value to cope with divergence 
of U at h = 0. 



D. Cluster size population kinetics 

Given the high ratio between the observed time scale 
of turbidity and the Brownian collision time of 4 % w/w 
silica particles, 3r?/(4fc s Tn ) < 10~ 4 sec pjj, the cluster 
growth under consideration is certainly RLA. 

In the most general scheme, both RLA and DLA clus- 
ter size population kinetics can be characterized by dy- 
namic scaling: X n — M f(n/M) where X n is the frac- 
tion of clusters of size n, and M(T) is an increasing clus- 
ter size characteristic of the system at a given time T 
[12, HH, where T — t/t agg is time scaled by the the initial 
monomer-monomer aggregation rate. The RLA cluster 
size dynamics is known to have a power law distribution 
described by the following three equations [H, [H, [24[ : 
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(5) 



7 



6 





ImM 


1 

lOmM 


1 1 
100 mM 


400mM _ 


















140mMNH 3 


buffer 






I 


I 





5 II i I i I i I i I i I 

0.5 1 1.5 2 2.5 



rj, e/nm' 



FIG. 1: pH control of a at different monovalent ionic 
strengths according to the Stern model of 20 nm silica. With 
4 % w/w silica, the buffering strength of 140 mM ammo- 
nia (Kb = 1.77 x 10 _5 M) is adequate to prevent significant 
change in pH with ionic strength. Details are explained in 
Appendix [A] 



M(T) = rVC 1 -*) (6) 

X n = C(T)(M(T)/n) 1+x (7) 

Eqn. ([S]) is the constraint of mass conservation. M repre- 
sents the number of particles in the biggest cluster of the 
entire system; it evolves in time according to Eqn. 
The third equation shows that X n depends on n in a 
power-law fashion with the characteristic exponent 1 + A. 
The scaling exponent A unites the evolution of both X n 
and the the size limit M. C(T) is chosen to satisfy 
mass conservation. Given a power law distribution of 
X n , the continuum integral 1 = nX n dn shows that 
C(T) = M~ 2 ; it can also be seen that M is proportional 
to the cluster size. This is summarized graphically in 
Fig. [2]for A = 0.5. Eqns. {7} and ([6]) can be derived from 
a scaling approach to the cluster aggregation process [23j . 

Simulation and experiment in the RLA regime show 
0.5 < A < 1.0 pi M, HI HE H3 and d f ~ 2.1 [H HI- 
The asymptotic limit of the scaling theory [24[ gives A = 
0.5. 
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III. MATERIALS AND METHODS 



X 




10 100 
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FIG. 2: A graphical representation of Eqns. (©-0 showing 
the self-similar nature of the time evolution of X n . In this 
example A = 0.5. M = T 2 (Eqn. (©) sets the upper limit of 
n for each instance of the distribution X n (T). 



Power law growth of M is seen experimentally [12], |27j 
and by simulations [2f| [2|| . There is a delay before the 
power law growth sets in, but the resulting M(T) may 
be time-scaled. 



E. Light scattered by a fractal cluster 

For an individual cluster of size n the total scattered 
light of wavelength A is 



S(k,R g )dk 



(8) 



where S(k, R g ) is scattering structure factor [2j|,|30j and 
the scattering vector k = (Att/X) sin(fl/2). For a cluster 
of fractal dimension d f I2al2al30|,l3l( 



S(k,R g 



{R g k) 2 ^-df/2 



3d//2 



(9) 



Thus the total scattering increases and the low angle scat- 
tering becomes greater as the cluster size R g becomes 
larger. 

The normalized scattered light r from a volume ele- 
ment of sample (the turbidity) would then be 



M 



B^Xr 



ir I S(k,R g {n))dk (10) 
feo 



where B is a constant [29[ . The spacial derivative of the 
light intensity dl/dx ~ —tI, thus the transmitted light 
follows a Beer's Law behavior, I ~ Jq exp(— tAx) for 
arbitrary Ax. 



In order to follow the aggregation of many samples 
of colloidal silica over a long time course, a carrousel 
was constructed to carry up to 32 standard 3 mL plastic 
sample cuvettes with a 1 cm optical path length. One 
cuvette served as the water reference. An incandescent 
bulb served at the measuring light source. Light was col- 
lected with a lens and focused into a pencil which passed 
through the rotating cuvettes. The transmitted light was 
detected with a photodiode; a slit in front of this pho- 
todetector made a 30 mrad angle of acceptance. A pi- 
coammeter recorded the photodiode current and was in- 
terfaced with a personal computer. The computer also 
controlled the rotation of the carrousel, allowing the light 
transmission to be recorded for a chosen time schedule. 
Temperature was controlled to be slightly above ambient, 
30° C. 20 nm colloidal silica was purchased from Alfa- 
Aesar (stock number 12727). PEG (Sigma 03557 and 
Alfa-Aesar A16151) was used as received from the ven- 
dors. Chemicals were of reagent grade. 

To initiate aggregation, 3M KC1 was introduced by 
slow addition (0.1 mL/min) with a syringe pump to the 
~ 3 mL volume of silica while being throughly mixed 
with a magnetically driven stirring propeller. A glass 
test tube of 1 cm inside diameter was used to contain 
the solution while mixing. The solution was then trans- 
fered to a plastic cuvette and capped with tape. 140 mM 
ammonia was used as the buffer for both the silica solu- 
tion and the 3M KC1 stock solution. Silica concentration 
was 4 % w/w. By investigating slow aggregation, prac- 
tical concerns about the exact initial conditions and the 
mixing process were minimized. The samples were never 
mechanically disturbed or shaken after the initial mixing. 
When PEG was part of the solution, it was added before 
the salt. 

A table-top demonstration of the acceleration of aggre- 
gation by PEG can be done by 1:1 v/v mixing of 4 % w/w 
silica with 20 % w/w 35 kD PEG. Both solutions have 
200 mM KC1 and 140 mM ammonia buffer. This silica 
suspension is undergoing aggregation, but very slowly. 
By layering a lighter PEG solution on top of a denser sil- 
ica solution, a whitish band of rapidly aggregating silica 
is observed at the interface. 

The diffusion constant, which is inversely proportional 
to viscosity, plays a role in the Smoluchowski model of 
aggregation. Viscosity data is available for 2 kD to 6 kD 
PEG in the literature [3^, but the range of concentra- 
tions measured do not extend into the semidilute Q 
regime c > c* . Because our measurements extended into 
this concentration range for 35 kD PEG, whose semidi- 
lute concentration starts at about 0.9 % w/w, the specific 
viscosities for PEG 4 kD and 35 kD PEG were measured. 
Viscosity measurements were performed with a 30 gauge 
stainless steel flow tube at 25°C. Water, 23 % w/w su- 
crose, and 46 % w/w sucrose served as calibration stan- 
dards. Our results are consistent with the published data 
f32| . There was no change in the quadratic trend of the 



specific viscosity for concentration of 35 kD PEG in the 
range to 4.5 % w/w. The relative viscosity for 35 kD 
can be fit with jj t = 1.0 + 0.49c + 0.124c 2 and 4 kD PEG 
with the equation r\ r = 1.0 + 0.166c where c is the % w/w 
concentration in the range to 4.5 %. 

The radially symmetric PB equation was numeri- 
cally integrated as previously described [33J. The so- 
lutions were found to agree with Eqn. (|A1|) for rele- 
vant values of pH and surface charge. The routines 
D01GAF and D01AJF from The Numerical Algorithms 
Group (http://www.nag.co.uk/, Oxford UK) were used 
to perform numerical integrations. 



KCl = 308 mM 
*o%peg 4 kD PEG 




IV. RESULTS 

A. Ionic strength and aggregation kinetics 

Fig. [3] presents a data set illustrating the effects of 
ionic strength on the aggregation kinetics of 20 nm silica. 
Each time course of the light transmission is empirically 
described with the function ao + aie _ (*/*°) 




KCl, mM 



FIG. 3: (a): Salt-induced aggregation followed by change in 
light transmission for 4 % w/w 20 nm silica. [KCl] = 260, 308, 
350, and 400 mM. Points are the measured light transmission, 
lines are a fit of the form ao+ai exp(— (t/to) 2,6 ). (b): The 1/to 
values of the four curves of (a) show a power law dependence 
on [KCl] with slope 6.5. 
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FIG. 4: The effect of 4 kD and 35 kD PEG on the aggregation 
kinetics with 308 mM ionic strength, 4 % w/w 20 nm silica, 
(a): 4 kD PEG [0, 0.87, 1.75, 2.6, and 3.5 % w/w]; (b): 35 kD 
PEG [0, 0.375, 0.747, 1.12, and 1.49 % w/w]. Points are the 
data, and lines are the parametric fit using the same stretched 
exponential function as in Fig. [3] 



We use the salt-induced aggregation time courses, 
shown in Fig. [3]Ja), as a means to evaluate our under- 
standing of the interaction U(r). A log representation 
of the dependence of l/t on ionic strength (Fig. [3{b)) 
demonstrates a power law with an exponent value of 6.5. 
Such a power-law relationship has been known for some 
time to exist over a limited range of salt concentration 
and is explained by the DLVO theory [101 ]. 

Eqn. (U) predicts the initial rate of dimer formation to 
be a function of n 2 . We do observe a quadratic depen- 
dence of 1/to on silica concentration (data not shown). 



B. Polymer-mediated interaction 

By setting the aggregation rate at a practical value 
with ionic strength, we can evaluate polymer-mediated 
effects. In Fig. UJa) and (b) we demonstrate the faster 
development of turbidity caused by 4 kD and 35 kD PEG, 
respectively. Notably, the stretched exponential power 
term a = 2.6 is preserved for all polymer concentra- 
tions. However, the final turbidity increases with increas- 
ing polymer concentration, whereas the final turbidity for 



6 



the salt-only data is almost constant. Clearly 35 kD PEG 
is more effective than 4 kD PEG in accelerating aggrega- 
tion. 



V. DISCUSSION 

A. Kinetics and the two-particle interaction 

As the integrand in Eqn. (J4| is of the form 
exp(J7(r)/fcsT), it is expected that the electrostatic in- 
teraction will have a strong effect on the aggregation rate. 
The dramatic effect of scaling the interaction by a con- 
stant is seen in Fig.^a). It is found empirically that an 
adjusted interaction 



(a) 



U x = 0.6 [£7, 



ele 



u, 



vdW 



(11) 



demonstrates a reasonable agreement with the 1/to val- 
ues of Fig.[3fb). These numerical integrations of Eqn. ((4]) 
use the interaction potential t/ eie + U v d w, with U e i e based 
on the values for a derived from Fig. [TJ A cutoff value 
c = 0.2 nm was used in Eqn. ((¥]) as U vc iw diverges at 
r = 2a. We use a Hammaker constant A = 1.66 fc^T, 
which is close to the commonly accepted Hammaker con- 
stant for silica particles in water of 2 ksT [Tlj |. 

The depletion interactions predicted by the PRISM in- 
tegral equation theory were added to U\, and the result- 
ing J compared with the experimental values of 1 /to from 
Fig. [TJ Calculations were also performed using the the 
RG theory according to Tuinier et aL[8|. Similar results 
were obtained (not shown). 

Both the RG and PRISM theories seem to overesti- 
mate the interaction. The experimental results may be 
matched by using 0.6 as a constant scaling factor. In 
summary 



U(r)=0.6[U v 
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FIG. 5: (a): The rate of dimer formation J is extremely 
sensitive to the interaction potential. By scaling Ui(r) by a 
factor 7 = 0.6, J matches the experimental values 1/to from 
Fig- El (b) : Calculation of J e / / using PRISM results. Data 
from Fig. [4] is included. The results of the calculation were 
scaled to match the data at zero polymer concentration. The 
PRISM interaction is from Ref. p], equations (8) and (12). 



J eff = J[U{r)]/ Vr (13) 

where J[U(r)] is numerically obtained from Eqn. Q. 
The effective rate of aggregation J e ff includes a factor 
of the relative viscosity, which is a function of the PEG 
concentration. In Fig. 0Jb) there is a rough agreement 
between the calculated aggregation rate and the mea- 
surement. 

Because there is no indication of Kramer inversion Q 
we believe that using macroscopic relative viscosity on 
the microscopic scale is proper. 



of integration k = 0, ki = 27r/4, df = 2.2, and A = 0.70, 
with the cluster size distribution of Eqns. (O, ([7|), and ([6]) 
yields the light transmission exp(— r), This calculation, 
shown in Fig. [SI closely resembles a stretched exponential 
form with a = 2.6 quite well. Various stretched exponen- 
tials may be obtained from different values of A and df 
(Table [TJ . We find the kinetics of light transmission is 
set primarily by A, not df. This can be understood by 
evaluation of the term J S(k,R g ) dk in Eqn. fTO")) . This 
integral is much more sensitive to R g than df, and the 
distribution of R g is set by A. 



B. Origin of the stretched exponential 

The stretched exponential decay of light transmission 
can be explained by combining the analysis of parts C 
and D. Numerical evaluation of r, Eqn. (|10p . with limits 



The calculations indicate zero light transmission as the 
aggregation proceeds, whereas the data, especially the 
salt-only data, show a limit in turbidity even after gela- 
tion is reached. This is certainly diffusion of photons in 
the medium [34| which we do not take into account. 

We have used three time scales here, to is the empir- 
ical fitting parameter for the measured light transmis- 
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FIG. 6: Light transmission calculation with A = 0.70 and 
df = 2.2 is able to reproduce the stretched exponential form 
of the data (exp(-(T/To) a )). In this case To = 5.16 and 
a — 2.6. The kinetics for other values of A and df are listed 
in Table U Light transmission is proportional to exp(— r). 

TABLE I: effect of df and A on fit parameters To and a for 
the light transmission fitting function exp _< - T//T °" ) . 





A = 0.65 


A = 0.70 


A = 0.75 


df = 2.1 


a = 2.21 
To = 6.72 


a = 2.55 
T = 5.25 


a = 3.01 
T = 4.07 


d f = 2.2 




a = 2.62 
T = 5.16 




df = 2.3 


a = 2.33 
To = 6.5 


a = 2.68 
To = 5.1 


a = 3.18 
To = 3.9 



sion; T is the dimensionless time of the dynamic cluster 
growth model; and then 1/J is the dimer formation time 
scale from the DLVO/Smoluchowski theory. Simulation 
[2! ] as well as numerical evaluation of the scaling model 
(Eqns. ©-©) indicates X 2 ~ 1/T, i.e., the time scale 
of dimer formation is basically T . So we can reasonably 
conclude T/5 = t/5t agg ~ t/to. The extreme sensitivity 
of J to the two-particle interaction potential (as seen in 
Fig. [5]) implies a rough accuracy in time scale is adequate 
to characterize the physics, so we can treat to and t agg 
as the same. 



centration (Fig. [4j. But remarkably the interplay of A 
and d f seems to preserve the time-scaling of the turbid- 
ity. According to Ball et al. [35], A is stabilized by the 
adjustment of df. To paraphrase: Imagine an increase 
in A, which leads to relatively more small clusters, as 
stated in Eqn. ([7]). These smaller clusters interpenetrate 
the larger clusters, which combine to form more compact 
objects of higher df. The resulting decrease in available 
surface area slows the growth, i.e., decreases A, which is 
consistent with Eqn. ([5]). 



D. Flavors of RLA 

Experimentally A = 0.5 has been found with very small 
colloid volume fractions, but with low colloidal stability 
set by either divalent ions or high concentration of NaCl 
[H,[l3|]. In these experiments, clusters seldom meet, but 
combine when they do meet. In our case, and in other 
experiments [13, HU, NaCl is about 100 mM, so clusters 
can meet many times before combining, and we find A = 
0.70. Two flavors of RLA are consistent with the analysis 
of Meakin and Family [26[ . 

It should be noted that 20 nm silica particles have a 
"hairy" nature presumably due to dangling poly(silicilic 
acid) chains on the surface [18j . Larger silica particles 
are seen to be " harder" and aggregates formed from such 
particles show a well-defined particle morphology [l8l ]. 
One might expect these smaller silica particles to form 
"large floppy clusters" *26j which have A = 0.70. 



VI. CONCLUSIONS 

To our knowledge this is the first study to combine the 
classical DLVO interpretation of colloidal stability with 
cluster-size dynamic scaling to further understanding of 
the AO depletion interaction. 

We are able to explain both the stretched exponen- 
tial nature of the light transmission data and the AO 
depletion interaction by qualitative adaption of current 
theories. 



C. Polymer effects and A 

For salt-induced aggregation (Fig. [3]) it appears the 
process is primarily characterized by to ■ Presumably the 
geometry of the aggregation (the cluster distribution and 
fractal dimension) at a given turbidity is the same for 
all ionic strengths; the ionic strength is simply setting 
the rate at which the process is played out. The nearly 
constant final turbidity implies the sample is proceeding 
to one common state regardless of salt concentration. 

The effect of PEG is more complex. Whereas salt- 
induced aggregation has only a small difference in tur- 
bidity at long time (Fig. [3}, PEG- induced aggregation 
shows increased final turbidity with increased PEG con- 
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APPENDIX A: THE STERN MODEL 

Here we explain in detail a model which explains the 
titration data of colloidal silica as typified by Fig. [TJ It 
combines the PB theory of electrolyte solutions, the Stern 
theory of surface capacitance, and the lattice model of the 
silica surface. 
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In a solution of monovalent salt of concentration no, 
the electrostatic potential surrounding a body with sur- 
face charge a can be found by integrating the PB equa- 
tion: 

V 2 4> = 2n o esinh(0e/fc B r)/e. 

with the boundary conditions a/eeo = d(f>(a)/dr and 
4>(oo) = 0. By solving this equation, we can find the 
surface charge as a function of surface potential <pd — 
4>(r = a), particle radius a, and Debye screening length 



m [21|: 

e 

2 

+ — tanh(e0 d /4fc B T)l 
na 

Because a monotonically increases with <f>d, we can also 
readily obtain (f>d(cr). 




FIG. 7: The Stern scheme of surface electrostatics. The 
(negative) surface charge a experiences potential (f>o. There 
is a potential drop a/Cs in going from the surface proper 
to the "diffuse" region where the counter-charge dwells. The 
potential experienced in the electrolyte medium is described 
strictly by <p(r = a) = (f> d and surface charge a. The Stern 
layer is merely characterized by the capacitance/area Cs- It 
can be of arbitrary thickness. 

Without specifying an exact physical mechanism, the 
Stern model proposes a capacitance Cs between the inner 
charged layer and the solution phase. The potential </>o 
of the surface charge is then 



<pd is the apparent surface potential experienced by the 
surrounding solution; 4>o is the potential of the charged 
SiO~ surface groups. 

We find it instructive to write the free energy F for one 
silica sphere in equilibrium with an ideal gas of protons 
in solution 1361: 



F — E e l e + Ebi n d 
pq=e(N — n) 



TS 

<}>o(q)dq + nu 
Nl 



(A3) 



■<r/Q 



s- 



(A2) 



9=0 

-Tk B log 

(Jy — n)\n\ 

+ k B T(M - n) log(((M - n)/V)/H ). 

where there are [N — n) SiO~ sites on the silica surface; 
n is the number of bound protons forming silanol (SiOH) 
sites. The successive terms are: the electrostatic charging 
energy; the binding energy for n protons, each forming 
a silanol group; the configuration entropy for n indistin- 
guishable protons distributed on N possible surface sites; 
and lastly, M — n protons in solution are treated as an 
ideal gas of volume V. This is simply Langmuir adsorb- 
tion with an electrostatic term. F is clearly a function 
of n, and by minimizing F(n) we find the surface charge 
(N - n)e. 

This Stern model, with free parameters density of sur- 
face sites N/4na 2 , proton binding energy uq, and sur- 
face capacitance Cs, is quite effective in explaining the 
pH titration data of colloidal silica particles for varying 
monovalent ionic strength [It], Ell]. The results shown 
in Fig. [1] represent the titration curves for several mono- 
valent ionic strengths obtained by numerically minimiz- 
ing F using the accepted values of 8 silane sites/nm 2 , 
u = pK k B T = 7.5 k B T, and C s = 2.9 F/m 2 [13, El- 

The energetics of the silica charged surface is typically 
dealt with in terms of chemical potential [2l[ ■ At equi- 
librium, the chemical potentials of a proton in the gas 
(solution) phase and a proton on the silica surface are 
identical, or equivalently dF/dn = 0. Also, the practical 
terms pH = — log 10 [H + ] and pK = uo/k B T are normally 
employed, hence the term "one-pK Stern model." 

In summary, for the given pH and ionic strength no, 
we can use the Stern model to find <j)d and in turn 
calculate the particle-particle electrostatic interaction 
U ele of Eqn. ©. 
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